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A flexible scheme for decomposing the vibrational density of states in terms of pair vibrational 
density of states is presented. This scheme provides the linkage between the site vibrational density 
of states and pair vibrational density of states so that vibrational modes, in particular localized 
modes, can be conveniently examined in terms of the correlation between the vibration at a given 
site and those at its neighboring sites. Furthermore, within the framework of a total energy vi- 
brational spectrum study, this scheme allows the analysis of vibrational modes in terms of their 
electronic origin. A case study of the vibrational dynamics of the relaxed Sig7 cluster is carried 
out to demonstrate the flexibility of the scheme in analyzing the properties of vibrational modes, 
particularly for complex systems with reduced or no symmetry. 

PACS numbers: 63.10.+a, 63.22.+m 



I. INTRODUCTION 



II. METHODS 



The dynamics of a system is described by the equation 
of motion 



In recent years, total-energy-vibrational-spectrum 
studies have been carried out for extended periodic sys- 
tems as well as complex systems with reduced symmetry 
because of the development of efficient and reliable meth- 
ods of calculating the total energy pi 0, 0- These stud- 
ies start with the calculation of the total energy using 
semi-empirical tight-binding (TB) approaches^ or first 
principles methods 0,0 such as those based on the den- 
sity functional theory (DFT). The force constant matrix 
between a pair of atoms are then calculated as the appro- 
priate second partial derivatives of the total energy with 
respect to the displacements of the atoms about their 
equilibrium positions. While it is indeed more satisfying 
to calculate the vibrational properties using force con- 
stant matrices determined from the total energy based 
on quantum mechanical calculations rather than from 
empirically-fitted classical potential, a key advantage of 
the total energy approaches over the approaches based 
on the classical potentials is that the former allows the 
possibility of linking the interesting dynamical features, 
in particular vibrational modes associated with surfaces, 
interfaces, defects, impuritites, etc., directly to the local 
electronic structure. However, this issue has not yet been 
systematically addressed. We present here a scheme that 
provides a framework through which dynamic features of 
a vibrating system can be related to the electronic struc- 
ture. We also demonstrate the application of this scheme 
using the vibration dynamics of the relaxed Sis7 cluster 
as the working example. We will show how all the in- 
teresting dynamical features of the relaxed Sis7 clusters 
can be traced to their electronic origins, including the 
local electronic structure such as bond charges, dangling 
bonds, etc. 
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where H ia . jP = <pi a ,jp/ y /m t m jl cf a 
fiia^p being the force constant matrix element between 
atom i of mass m, along the direction a (x,y, or z) 
and atom j of mass nij along the direction (3, and uf a 
the eigen-displacement of the A-mode (with the eigen- 
frequency u>\) at site i along the a-direction. Tradi- 
tionally, the vibrational density of states (VDOS) pro- 
jected on sites (individual atoms) is the tool to analyze 
the nature and properties of a defect or impurity mode. 
This is usually accomplished by introducing the sum rule 
J2i a ( c ia) 2 = 1 m t° the defining equation of the VDOS 
such that 

pM = £ 5(u - u x ) = Y,(ct) 2 5(u - w x ) (2) 
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In this way, 



where 



,H = X)(«fa) a <y(w-«A) 
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is the local VDOS at site i along the a polarization. How- 
ever, since the force field characterizing a particular vi- 
brational mode is determined by the interaction of the 
atom at a given site with its neighbors, it is more illu- 
minating if, in addition, the VDOS projected on pairs of 
atoms is used in combination with the site VDOS to ana- 
lyze the mode. The interactions of a certain atom with its 
neighbors are described by the force constant matrices. 
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Thus a sum rule that involves the force constant matrices 
would be an appropriate vehicle to construct the VDOS 
projected onto pairs of atoms. Such a scheme can be set 
up as follows. 

Multiplying the equation of motion with c\ a (Eq. (1)), 
we have |1 0| 



(cL) 2 = E 
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In arriving at Eq. (5), the relation <j>i a> ip = ~ J2j=^i 4>ia,j0 
has been used. Equation (5) can serve two useful func- 
tions. First, the summation of Eq. (5) over both i and a 
leads to 

( u ia u jl3 ~ u ioi U ip)4'ioi,30 _ \^ i „X \2 



(<&r = i (6) 



This is precisely the sum rule needed to project the 
VDOS onto pairs of atoms i and j as Eq. (6) also in- 
volves the force constant matrix elements 4>i a ,j/3- Thus, 
when Eq. (6) is inserted into the defining equation for 
the VDOS (Eq. (2)), we obtain 
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where Pi a ,j p(w), the pair VDOS of atom i along the in- 
direction and atom j along the /3-direction linked by the 
force constant matrix element <fri a ,jp, is given by 

puxMu) = 2^ 2 — <n w - w a) (8) 

X W A 

The second useful function of Eq. (5) is that it provides 
the framework to link the local VDOS at a particular 
site to the interactions with its neighbors. In this way, 
the vibrational mode localized in the vicinity of a given 
site can be correlated with the vibrational motions of its 
interacting neighbors. Thus, from Eqs. (4) and (8), we 
have 



Equation (5) can also be rewritten as 
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1 then leads to 



PH = 2J Pi<*,30( U) ) 
ia,j0 
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with 

, , ^ <t>ia,30 u ia u ja - <t>i a ,ia( u D 2S i3 § <xP £ , , 
Pia,j0{U) = 2^ 2 _ tt 0{CJ-CJX) 
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Equations (10) and (12) also yield 



Pia(to) = 2J Pia,j0(w) 
30 



(13) 



Equation (13) allows the decomposition of the site VDOS 
in terms of on-site as well as off-site contributions. 

While the scheme based on Eq. (5) provides the ve- 
hicle to decompose the VDOS in terms of contributions 
from pairs of atoms directly interacting with each other 
via 4>i a ,jf3, one may also devise a scheme to project the 
VDOS onto pairs of atoms interacting directly as well as 
those not interacting directly. Such a scheme is useful for 
examining collective (extended) as well as more localized 
vibrational modes. To develop such a scheme, we note 
that 
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Multiplying Eq. (14) by c; la and summing over i 7 a and 
(3 lead to 

E c »4? = E( c -) 2 + E ^0 - (4j%<m = 
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Equation (16) then leads to 



PH = E Pi<*jP( U ) 
ia,j0 



(17) 



with 



Pia,j0(u) = - E{ C «* C J/3 " (4J 2 <%<W*M( W ~ W A) (18) 



According to Eq.(16), we can also write 

Pia(uj) = E Piajp(v) 
30 
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This equation then provides the link of the site VDOS to 
the pair VDOS Pi a ,jp- Similarly J2j c jp — leads to 

E c »4* = E( c -) 2 + E = (20) 



or 



p(uj) = E Piajafa) 



(21) 
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with 

ftaJaM = - <&<$J<V - ^) ( 22 ) 
A 

Equation (20) also leads to 

Pia(w) = y~]pia,ja(w) (23) 
3 

Equations (5) through (23) provide some examples of 
projecting the VDOS onto pairs of atom and relating 
site VDOS to pair VDOS for the purpose of studying the 
correlation between the vibrations at different localities 
for localized as well as more extended vibrational modes. 

Within the framework of a total-energy-vibrational- 
spectrum approach, the bonding configurations of atoms 
in the system, such as the bond charges, the nature of 
the bonds, and the number of bonds associated with each 
atom, are already determined [(j. Therefore, it is straight- 
forward to relate a particular vibrational mode to its elec- 
tronic origin via either the site VDOS or the pair VDOS. 
For example, the VDOS can be partitioned according to 
the contributions from different ranges of the total bond 
charge % = (1/2) JV cr^- associated with the atom at site 
i in the system where <7y is the bond charge between the 
pair of atoms i and j . This can be accomplished by sum- 
ming the site VDOS pi = pi a over those sites i where 
qi falls within the range [a, b], i.e., 

p(u; a < qi < b) = PiM ( 24 ) 

Since the bond charge is closely related to the strength of 
the bond, such a partitioning can shed light on whether 
a particular mode is associated with the strengthening or 
softening of the bonds. Similarly, the VDOS can also be 
grouped in terms of contributions from atoms possessing 
a certain number of bonds, namely 

p(w;m = k) = Pi(uj) (25) 

rii—k 

where rij denotes the number of bonds possessed by the 
particle i. This presentation of the vibrational spectrum 
provides the means to identify the modes associated with 
dangling bonds and those associated with distorted bond- 
ing configurations. To examine how an individual bond 
affects a particular mode, one may group the pair VDOSs 
Pij = YlaB Piajp according to the range of the bond 
charge ay such that 

p(oj; a < aij < b) = 2J PijM ( 26 ) 

a<(Tij<b 

The characteristics of the bonds is usually described 
by molecular orbitals such as r = ssa, spa, ppa, or ppn. 
To determine the nature of the bonds responsible for the 



vibrational modes, one may classify the contribution to 
a particular mode by a certain bonding type r using 

P {u;t)=Y j —pM (27) 
rrl. <J%j 

where 

(28) 

r 

The VDOS can also be decomposed in terms of quan- 
tities such as the distance, the bond charge, etc. via 
Pia,j/3, or Pia, j f3 or p~i a ,j a , following similar procedure as 
described in Eq. (26). Thus, a very flexible framework 
that allows the convenient decomposition of the VDOS 
to pair VDOS via p ia jp, or p Mj73 , or p ia ,jp, or p ia , ja , 
and the linking of the site VDOS to these pair VDOS (see 
Eqs. (9), (13), (19), and (23)) can be set up. This frame- 
work can then be profitably used to analyze interesting 
dynamical features of vibrations, in particular how these 
features are correlated with the local electronic structure, 
of complex systems with reduced or no symmetry. 

III. RESULTS AND DISCUSSION 

To demonstrate the application of the scheme outlined 
in Eqs. (5) to (27) to examine the properties of the vibra- 
tional modes and to trace these modes to their electronic 
origins, we use the dynamics of a relaxed Sis7 as the work- 
ing example. The initial configuration of the cluster is an 
87-atom tetrahedral network with an atom at the central 
site. The x-, y-, and z-axis are along the cube edges of 
the cube-unit cell with the central site at a corner of the 
cell. The equilibrium configuration of this cluster (see 
Fig. 1) was first determined by the MD scheme using 
the non-orthogonal tight-binding Hamiltonian developed 
by Mennon and Subbaswamy0. The force constant ma- 
trices of the relaxed cluster were then numerically calcu- 
lated from the total energy. 

In Fig. 2, the VDOS of the relaxed cluster, together 
with contributions to the VDOS from those sites with 
the total bond charge in a certain range (see Eq. (24)), 
is shown. We use the bonding configuration of the bulk 
Si as the reference to choose the ranges of the total bond 
charge because Si-based systems favor the tetrahedral sp 3 
bonding configuration. In bulk Si, each Si atom con- 
tributes 0.25e to one of the four equivalent bonds, giving 
rise to a total bond charge of le for every atom. Thus, we 
have selected the following four ranges, namely qi > le, 
le > qi > 0.8e, 0.8e > qi > 0.6e, and 0.6e > qi, to exam- 
ine the effect of relaxation on the bonding configuration, 
and hence on the vibrational spectrum of Sis7 . From Fig. 
2, it can be seen that the majority of the contributions to 
the VDOS originate from those sites with le > qi > 0.8e 
or 0.8e > g, > 0.6e, indicating that the relaxation has 
probably caused a general softening of the bonds and 
thus resulting in an overall shift of the VDOS towards 
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lower frequencies compared to the bulk VDOS. Figure 2 
also shows that the mode with the highest frequency is 
a defect mode originating from the range > le, thus a 
strengthened bonding configuration, while the mode with 
the lowest frequency originates from the range 0.6e > 
a softened bonding configuration or a configuration with 
dangling bonds. In Fig. 3, the contribution to the VDOS 
is decomposed according to the number of bonds associ- 
ated with each atom (see Eq. (25)). The existence of a 
bond between a pair of atoms is established by the cri- 
terion cry > 0.04e rather than relying on the distance^]- 
Figure. 3 shows that the majority of the contributions 
still originate from atoms with 4 bonds, with a signif- 
icant minority contribution from atoms with 3 bonds. 
The defect mode at the highest frequency can be seen as 
associated with an atom with 4 bonds. Together with the 
evidence in Fig. 2, it can be concluded that this mode re- 
sults from the strengthened bonding configuration. The 
lowest (non-zero) frequency mode, on the other hand, is 
seen to originate from an atom with 2 bonds. The ev- 
idences provided by Figs. 2 and 3 thus indicate that 
this mode is a dangling-bond mode, but with remain- 
ing bonds strengthened. To further explore this issue, 
we present in Fig. 4 the decomposition of the VDOS in 
terms of the ranges of bond charge cr^ (see Eq. (26)). 
ft can be seen that indeed both the highest- and the 
lowest-frequency mode are associated with <7jj > 0.5e, 
indicating that the bonds involved in these vibrations are 
strengthened bonds (compared to the bulk bond) . Figure 
4 also shows that the majority of the contribution to the 
VDOS come from the range 0.4e < cr^ < 0.5e, the range 
with bond charge not too much deviated from the bulk 
bonding configuration, thus reinforcing the observation 
provided by Figs. 2 and 3. Another interesting feature 
in Fig. 4 is the negative contribution in the low frequency 
VDOS from the range ov,- < 1 0~ 3 e. The negative contri- 
bution is an indication that for contributions to a certain 
vibrational mode by pairs of atoms i and j interacting 
via tpia.jfs, some of those with no chemical bond formed 
between them may have the force constant too "soft" to 
support the vibration at that frequency. In Fig. 5, we 
show the decomposition of the VDOS in terms of bond- 
ing types (see Eq. (27)). It can be clearly seen that 
the main contribution to all the vibrational modes comes 
from ppcr-bonds, with a substantial minority contribu- 
tion from sptr-bonds. The directional characteristics of 
the pper-bonds indicates the dominance of the stretching 
vibrations as is expected for a Si-based structure. 

Figure 6 presents the decomposition of the VDOS ac- 
cording to the ranges of the distance between atoms via 
Pia.jp (see Eqs. (18)). Since Pia.jp only involves cf a and 
c^g, it provides an useful vehicle to analyze contributions 
to a certain vibrational mode from the on-site environ- 
ment, pair of atoms directly interacting (via <fii a ,jp) with 
each other, and pairs of atoms not directly interacting. 
The three ranges of distances chosen for the decomposi- 
tion in Fig. 6 correspond roughly to 1st nearest-neighbor 
distances (r^ < 3. 10 A), 2nd nearest-neighbor distances 



(3. 10 A < rij < 4. 15 A), and 3rd nearest-neighbor dis- 
tances and beyond (r^ > 4.15A) respectively. Figure 6 
shows that there are contributions to the VDOS from 
all three ranges to most of the low- to intermediate- 
frequency modes and the contributions are consistent 
with the number of pairs of atoms in each range. This 
is an indication that these modes are extended modes. 
On the other hand, the dominant contribution to the 
high-frequency modes mainly originates from the range 
< 3.10A. Considering the fact that for any given 
atom in the cluster, there are many more neighbors at 
farther distances than close distances, the dominance of 
the contribution to VDOS from < 3. 10 A for high- 
frequency modes suggests that these modes are more lo- 
calized. For example, the highest-frequency mode ap- 
pears to be a localized mode as contributions to this 
mode from second nearest neighbors and beyond are out- 
weighted by the contribution from the nearest neighbors, 
even though there are many more farther neighbors than 
nearest neighbors for every atom in the cluster. The 
inset of Fig. 6 shows in more detail the contributions 
to the non-zero lowest-frequency mode. The dominant 
contribution to this mode is apparently from the on-site 
environment. To gain more insight to these two inter- 
esting modes, one may first plot J2 a ( c ia) 2 f° r the mode 
lj\ as a function of site index i to determine where the 
mode is localized and then decompose the site VDOS for 
that site according to distances (1st-, 2nd-, etc., based 
on Eqs. (8) and (9)) to shed light on the correlation be- 
tween the vibration of the site and those of its neighbors. 
The upper inset of Fig. 7 is the plot of J2 a ( c ia) 2 f° r the 
highest-frequency mode ( 108 THz) as a function of the 
site index. The mode is obviously localized in the imme- 
diate neighborhood of site 1 (the central site). The main 
panel of Fig. 7 shows the site VDOS p\ — ^ a Pia{u) 
decomposed in terms of 1st-, 2nd-, etc., nearest-neighbor 
distance using Eqs. (8) and (9). It can be seen that 
the contribution to this mode is almost entirely from the 
coupling of the central site (site 1) with its 4 nearest 
neighbors. In the lower inset of Fig. 7, the nature of 
the coupling is examined by decomposing p\ in terms of 
Y^jjti Pia,j/3- It is seen that the majority of the contribu- 
tion are from xx-, yy-, and zz-nearest neighbor coupling. 
It is also seen that the mode is actually a composite of 
three almost degenerate modes (107.2, 108.0, and 107.5 
THz), an indication of the distortion of the immediate 
environment of the central site from the perfect tetrahe- 
dral configuration. 

The analysis of the non-zero lowest-frequency mode is 
shown in Fig. 8. The upper inset shows that the mode is 
an on-site mode at site 82 located on the surface of the 
cluster (see Fig. 1). The site VDOS p$2, decomposed 
according to Eq. (13), confirms this picture. The lower 
inset of Fig. 8 is the decomposition of p$2 according to 
Eq. (13) so that the on-site coupling can be examined. 
It shows that two almost equivalent on-site xy- and yx- 
coupling are the main contributors to this mode, indicat- 
ing that this is mainly a wagging mode in the xy plane 
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anchored by the two strengthened bonds of the atom at 
site 82 (a site with two dangling bonds, see Fig. 1). 

The relaxed Sis7 cluster is a system with no symme- 
try The demonstration discussed in this section shows 
the flexibility of the methodologies developed in Sec. II 
to carry out detailed analysis of the dynamics of such 
a system and to trace interesting dynamical features to 
their electronic origin. Hence the framework described in 
Sec. II is expected to be extremely useful for the study 



and the analysis of vibrational dynamics of complex sys- 
tems with low or no symmetry 
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FIGURES 

FIG. 1 The equilibrium configuration of the relaxed Sis7 
cluster. The z-axis is the vertical axis and the x-axis is 
the horizontal axis. Atom 82 (see the result and discus- 
sion in the text) is highlighted by the dark ball. 

FIG. 2 The decomposition of the VDOS in terms of the 
range of the bond charge qi associated with the atoms in 
the cluster. 

FIG. 3 The decomposition of the VDOS in terms of the 
number of bonds associated with the atoms in the cluster. 

FIG. 4 The decomposition of the VDOS in terms of the 
bond charge cry between pairs of atoms in the cluster. 



FIG. 5 The decomposition of the VDOS according to the 
bond types. 



FIG. 6 The decomposition of the VDOS in terms of the 
distance between pairs of atoms according to Piajp (Eq. 
(17)). The inset is the enlargement of the curves in the 
immediate neighborhood of the mode with the lowest 
non-zero frequency. 

Fig. 7 Upper inset shows the plot of ^ Q (c^ Q ) 2 vs. the site 
index i for the highest frequency mode ( 108 THz). The 
main panel shows the site VDOS for site 1 (the central 
site), pi, decomposed according to 1st nearest neighbor, 
2nd nearest neighbor, etc., distance using Eqs. (8) and 
(9). The lower inset shows the decomposition of pi in 
terms of Piajp m the immediate neighborhood of 

the mode with the highest frequency. 



Fig. 8 Upper inset shows the plot of J2 a ( c ia) 2 vs - the 
site index for the lowest non-zero frequency mode. It can 
be seen that this mode is confined at the site 82 (a surface 
site, see Fig. 1). The main panel shows the site VDOS 
PS2 and its decomposition according to Eq. (13). The 
lower inset shows the decomposition of p 82 in terms of 
J2j Ps2a,j/3 in the immediate neighborhood of the mode 
with the lowest non-zero frequency. 
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